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We consider the time evolution of entanglement in a finite two dimensional transverse Ising model. 
The model consists of a set of 7 localized spin-i particles in a two dimensional triangular lattice 
. . , coupled through nearest neighbor exchange interaction in presence of an external time dependent 

■ magnetic field. The magnetic field is applied in difi'erent function forms: step, exponential, hyper- 

' bolic and periodic. We found that the time evolution of the entanglement shows an ergodic behavior 

under the effect of the time dependent magnetic fields. Also we found that while the step magnetic 
(N : field causes great disturbance to the system creating rabid oscillations, the system shows great 

f~| ' controllability under the effect of the other magnetic fields where the entanglement profile follows 

^ closely the shape of the applied field even with the same frequency for periodic fields. This follow 

. up trend breaks down as the strength of the field, the transition constant for exponential and hy- 

' perbolic, or frequency for periodic field increase leading to rapid oscillations. We observed that the 

I entanglement is very sensitive to the initial value of the applied periodic field, the smaller the initial 

, value the less distorted is the entanglement profile. Furthermore, the effect of thermal fiuctuations 

I I ' is very devastating to the entanglement which decays very rapidly as the temperature increases. 

' Interestingly, although large value of the magnetic field strength may yield small entanglement, it 

was found to be more persistent against thermal fiuctuations than the small field strengths. 
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C ' PACS numbers: 03.67.Mn, 03.65.Ud, 75.10.Jm 
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I. INTRODUCTION 

Quantum entanglement lies in the heart of quantum theory and is of fundamental role in modern physics [J]. 
Entanglement is a nonlocal correlation between two (or more) quantum systems such that the description of their states 
QQ ] has to be done with reference to each other even if they are spatially well separated. Understanding and quantifying 
, entanglement may provide an answer for many questions regarding the behavior of complex quantum systems [H. 
l/^ ' Particularly, entanglement is considered as the physical property responsible for the long-range quantum correlations 
' , accompanying a quantum phase transition in many-body systems at zero temperature (3|-l7|. Particular fields where 
■ entanglement plays a crucial role are quantum teleportation, quantum cryptography and quantum computing, where 
T-H [ it is considered as the physical basis for manipulating linear superpositions of the quantum states to implement the 
^-H . different proposed quantum computing algorithms [8[ Q . Different physical systems have been proposed as reliable 
^ ' candidates for the future technology of quantum computing and quantum information processing jlCll4l7| . The main 
task in each one of these systems is to specify certain quantum degree of freedom as the qubit and find a controllable 
coupling mechanism to form an entanglement among these qubits to perform efficient quantum computing processes. 

Multiparticle systems are of central interest in the field of quantum information, in particular, quantification of 
the entanglement contained in their quantum states. However, quantum states and entanglement are very fragile due 
to the induced decoherence caused by the inevitable coupling to the environment. Decoherence is considered as one 
of the main obstacles toward realizing an effective quantum computing system [18l] . The main effect of decoherence 
is to randomize the relative phases of the possible states of the considered system. Quantum error correction [l^ 
and decoherence free subspace j20l . [2l| have been proposed to protect the quantum property during the computation 
process. Still offering a potentially ideal protection against environmentally induced decoherence is difficult. In NMR 
quantum computers, a series of magnetic pulses were applied to a selected nucleus of a molecule to implement quantum 
gates [2^ . Moreover, a spin-pair entanglement is a reasonable measure for decoherence between the considered two- 
spin system and the environmental spins. The coupling between the system and its environment leads to decoherence 
in the system and vanishing of entanglement between the two spins. Evaluating the entanglement remaining in the 
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considered system helps us to understand the behavior of the decoherence between the considered two spins and their 
environment |23i] . 

In previous works, the evolution of entanglement in a one-dimensional spin system in presence of different forms of 
external magnetic fields, modeled by the XY Hamiltonian, was studied [23, ISa]- It was found that the entanglement 
can be localized between nearest-neighbor qubits for certain values of the external time-dependent magnetic fields. 
Moreover, as known for magnetization of this model, the entanglement showed nonergodic behavior, i.e. it does 
not approach its equilibrium value at the infinite time limit. Also, the same system was investigated considering 
a time-dependent exchange coupling between neighboring spins [26j . It was found that the asymptotic behavior of 
entanglement at the infinite time limit is very sensitive to the initial values of the coupling and magnetic field and 
for particular choices we may create finite asymptotic entanglement regardless of the final values of coupling and 
magnetic field. The quantum effects in the Ising model case showed persistence in the vicinity of both its critical 
phase transition point and zero temperature as it evolves in time. 

The study of quantum entanglement in two-dimensional systems possesses a number of extra problems compared 
with systems of one dimension. The particular one is the lack of exact solutions. The existence of exact solutions has 
contributed enormously to the understanding of the entanglement for ID systems [13, HB-Hl] ■ Studies can be carried 
out on interesting but complicated properties, applied to infinitely large system, and so forth use finite scaling method 
to eliminate the size effects, etc. Some approximation methods, like Density matrix renormalization group (DMRG), 
are also only workable in one dimension p2l . [29l - [3l| . So when we carry out this two-dimensional study, no methods 
can be inherited from previous researches. They heavily rely on numerical calculations, resulting in severe limitations 
on the system size and properties. For example, dynamics of the system is a computational-costing property. We 
have to think of a way to improve the effectiveness of computation in order to increase the size of research objects; 
mean while dig the physics in the observable systems. It may show the general physics, or tell us the direction of 
less resource-costing large scale calculations. In a previous work we have studied the entanglement in a 19-site two- 
dimensional transverse Ising model at zero temperature [33 |. The spin-i particles are coupled through an exchange 
interaction J and subject to an external time-independent magnetic field h. We demonstrated that for such a class 
of systems the entanglement can be tuned by varying the parameter X = h/ J and also by introducing impurities into 
the system. The system showed a quantum phase transition at a specific critical value of the parameter Ac. 

In this paper, we consider the dynamics of entanglement in a two-dimensional spin system, where spins are coupled 
through an exchange interaction and subject to an external time-dependent magnetic field. Four forms of time- 
dependent magnetic field are considered: step, exponential, hyperbolic and periodic. To tackle down the problem, we 
introduce two calculation methods: step by step time-evolution matrix transformation and step by step projection. 
We compare them side by side, in short, besides the exactly same results, step by step projection method turned 
out to be twenty times faster than the matrix transformation. One section of the paper will expound the scalability 
of this method, where it sheds light on studies of larger systems. The finite temperature effect is also considered 
to simulate more realistic systems. We show that the system entanglement behaves in an ergodic way in contrary 
to the one-dimensional Ising system. The system shows great controllability under all forms of external magnetic 
field except the step function one which creates rapidly oscillating entanglement. This controllability is shown to be 
breakable as the different magnetic field parameters increase. Also it will be shown that the mixing of even a few 
excited states by small thermal fluctuation is devastating to the entanglement of the ground state of the system. 
The critical temperature at which the concurrence vanishes depends significantly on the value of the magnetic field 
strength, smaller value yield smaller entanglement but higher critical temperature. 

This paper is organized as follows. In the next section we present our model and discuss the two different approaches 
to evaluate the entanglement. In sec. Ill we present and discuss our results for the entanglement of the system under 
the effect of the different magnetic fields forms at zero temperature. The dynamics of thermal entanglement is 
considered in sec. IV. In sec. V we explain part of our key results in the light of Fermi's golden rule and adiabatic 
approximation. The extension of our work to larger size spin systems is discussed in sec. V. We conclude in Sec. VI 
and discuss future directions. 



II. SOLUTION OF THE TIME-DEPENDENT TWO-DIMENSIONAL ISING MODEL 

A. Model 

We consider a set of 7 localized spin-^ particles in a two dimensional triangular lattice coupled through exchange 
interaction J and subject to an external time-dependent magnetic field of strength h{t). The Hamiltonian for such a 
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system is given by 

H = -Y.J.,,ata^^^h{t)Y,al, (1) 

<i;j> i 

where < ?,j > is a pair of nearest-neighbors sites on the lattice, J^j = J for aU sites. For this model it is convenient 
to define a dimensionless coupling constant A = h/ J . We apply different forms of the magnetic field as a function of 
time: step function, exponential, hyperbolic and periodic. 



B. The evolution operator 

According to quantum mechanics the transformation of \ipi{to)), the state vector at the initial instant to, hito \^i{t)), 
the state vector at an arbitrary instant, is linear [33|. Therefore there exists a linear operator U{t, to) such that: 

^uit,to)\Mto))- (2) 

This is, by definition, the evolution operator of the system. Substituting Eq. ([2]) into the Schorodinger equation, we 
obtain: 

ih-^U{t,to)mto)) = H(t)U{t,to)mo)), (3) 

which means 

ih^U{t,to) = H{t)U{t,to). (4) 

Further taking the initial condition 

U{to,to)=l, (5) 
the evolution operator can be condensed into a single integral equation: 

i /■* 

U{t,to)=I-- H{t')U{t',to)dt'. (6) 

When the operator H does not depend on time, Eq. ^ can easily be integrated and at finally gives out 

C/(t,io) = e-'^(*-*°)/^ (7) 



C. Step by step time-evolution matrix transformation 



To unveil the behavior of concurrence at time t, we need to find the density matrix of the system at that moment, 
which can be obtained from 

Pit) = U{t)p{Q)U\t). (8) 
Although Eq. ([ni gives a beautiful expression for the evolution operator, in reality U is hard to be obtained because 
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FIG. 1: The external magnetic field in a step function form h{t) = a @ t < to, h{t) = b @ t > tg. 

of the integration involved. In order to overcome this obstacle, let us first consider the simplest time-dependent 
magnetic field: a step function of the form (Fig. [T]) 



h{t) = a + {b - a)d{t - to) , 



(9) 
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where 9(t — to) is the usual mathematical step function defined by 

Oit-to) = [[ (10) 

At to and before, the system is time-independent since Ha = H(t < to) = — J2<i j> '^f'^j ~ aYnaf. Therefore we 
are capable of evaluating its ground state and density matrix at straightforwardly. For the interval to to t, the 
Hamiltonian Hb = H{t > to) — — J2<i j> '^f'^j ^ bT,iaf does not depend on time either, so Eq. ([7]) enables us to 
write out 

t/(i,to) = e''^(*>*'')(*-*°)/^, (11) 

and therefore 

p{t) = U{t,to)p{to)U^{t,to). (12) 

Starting from here, it is not hard to think of breaking an arbitrary magnetic function into small time intervals, and 
treating every neighboring intervals as a step function. Comparing the two graphs in Fig. [2l the method has just 
turned a ski sliding into a mountain climbing. Assuming each time interval is At, setting h = 1 then 
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FIG. 2: (Color online) Divide an arbitrary magnetic field function into small time intervals. Every time step is At. Treat the field within 
the interval as a constant. In all, we turn a smooth function into a collection of step functions, which makes the calculation of dynamics 
possible. 

U{t„to)\4'o) = U{t^,t,-l)U{t^-l,t^-2)■■■U{tl,to)\^o), (13) 

i 

U{t^,to) = X{exp[-iH{tk)M], (14) 
fe=i 

U{t„to) = exp[~iH{ti)At\U{t,_^-to). (15) 

Here we avoided integration, instead we have chain multiplications which can be easily realized as loops in computa- 
tional calculations. This is a common numerical technique; desired precisions can be achieved via proper time step 
length adjustment. 



D. Step by step projection 

Step by step matrix transformation method successfully breaks down the integration, but still involves matrix 
exponential, which is numerically resource costing. We propose a projection method to accelerate the calculations. 
Let us look at the step magnetic field again Fig. [T] For Ha, after enough long time, the system at zero temperature 
is in the ground state |(/)) with energy, say, e. We want to ask how will this state evolves after the magnetic field is 
turned to the value b? Assuming the new Hamiltonian Hi, has N eigenpairs Ei and \il>i). The original state 10) can 
be expanded in the basis {iV'i}}- 

\4>) = cilVi) + C2\^2) + - + catIV-at), (16) 

where 



Ci = {tpi\(t)). 



(17) 
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When H is independent of time between t and to then we can write 

U{t, to) IVmo) - e-^^(*>*'')(*~*°^/''|^Mo) - e-'^'(*-*'')/''|^,.t„), (18) 

Now the exponent in the evolution operator is a number no longer a matrix. The ground state will evolve with time 
as 

N 

= ^ c, I V',)e-^^' (19) 

4=1 

and the pure state density matrix becomes 

pit)^\m{m- (20) 

Again any complicated function can be treated as a collection of step functions. When the state evolves to the next 
step just repeat the procedure to get the following results. Our test shows, for the same magnetic field both methods 
give the same results, but projection is much (about 20 times faster) than matrix transformation. This is a great 
advantage when the system size increases. But this is not the end of the problem. The summation is over all the 
eigenstates. Extending one layer out to 19 sites, fully diagonalizing the 2^^ by 2^^ Hamiltonian and summing over all 
of them in every time step is breath taking. In a later section, we will show how to further improve the method to 
shrink the amount of calculation which will pave the way towards larger systems. 



E. Entanglement of formation 



We confine our interest to the entanglement of two s pins , at any position i and j [3J]. We adopt the entanglement 
of formation, a well known measure of entanglement [35l |. to quantify our entanglement [sgI ]. All the information 
needed in this case, at any moment t, is contained in the reduced density matrix Pi.j{t) which can be obtained form 
the entire system density matrix by integrating out all the spins states except i and j. Wootters [35| has shown that, 
for a pair of binary qubits, the concurrence C, which goes from to 1, can be taken as a measure of entanglement. 
The concurrence between sites i and j is defined as 

C{p) = max{Q, ei - 62 - 63 - £4}, (21) 

where the e^'s are the eigenvalues of the Hermitian matrix R = yf^Jpp^/p with p — [a^ ® a'^)p*{ay (g) cr^) and cr^ is 
the Pauli matrix of the spin in y direction. For a pair of qubits the entanglement can be written as, 

E{p) = e{C{p)), (22) 

where e is a function of the "concurrence" C 

<C) = h (^^4^) ' ^2^^ 

where h is the binary entropy function 

h{x) = — a;log2 a; — (1 — x)log2{\ — x). (24) 

In this case, the entanglement of formation is given in terms of another entanglement measure, the concurrence C. The 
matrix elements of the reduced density matrix needed for calculating the concurrence can be obtained numerically 
using one of the two methods developed above. 



III. DYNAMICS OF THE SPIN SYSTEM IN A TIME-DEPENDENT MAGNETIC FIELD 



A. Step magnetic field 



At first, we study the dynamic response of the 7-site spin system to the external step magnetic field (Fig. [T]), the 
simplest form of time-dependent function, which is given by 

h{t) = a + {h - a)e{t ~ to) , (25) 



6 



In a previous work we have showed that this two dimensional spin system has a critical point at A = 1.64 and a 
maximum reachable concurrence at A = 2.61 (32| . Accordingly, we design here all kind of steps: big, small, jump, 
drop, happening before, after or cross the "critical point" (and "maximum point"), trying to identify what will affect 
the system behavior and how. 




FIG. 3: (Color online) Dynamics of C(l,4) in the 7-site system when the step magnetic field is changed from a = 0.5 to fe = 1.5 (before 
the "critical point" h = 1.64), from a = 2 to fe = 3 (after) and from a = 0.5 to fe = 3 (big step cross the "critical point"). 

Figure [3] displays the dynamics of the pairwise entanglement between sites 1 and 4 C(l,4) when the magnetic field 
is changed from a = 0.5 to 6 = 1.5 (before the "critical point" h = 1-64), from a = 2 to b = 3 (after) and from 
a — 0.5 to b — 3 (big step cross the "critical point"). Thick and big oscillations appear in every graph. Dynamics of 
C(l,4) under similar design (Figure[5|), but before, after and cross the "maximum point" h = 2.61, shows resembling 
oscillations. In these graphs, we have plotted the concurrence corresponding to a constant magnetic field h = a (red 
line), h — b (yellow line) and the average value of the oscillating concurrence (green line). In the same manner, the 
dynamics of the concurrence C(l,2) is explored in Figs. S] and IH] which shows a very similar behavior to the C(l,4) 
case. Examining the dynamics of the next nearest neighbor concurrences C(l,5) and C(l,7) in the same way as we 
did with C(l,2) and C(l,4), we observed a very similar behavior with much much smaller value of the concurrence as 
expected for next nearest neighbors. 

Figure [7] intends to check behaviors of the system going from one magnetic field value to a smaller one and the 
reverse process. Both of them show the same characteristic oscillations, although they are different. Because the 
projection method requires mapping the ground state h\ h = a into all the eigenstates of the system in /i = 6; while 
we reverse the process, mapping becomes from the ground state \\\ h — b to all the eigenstates of the system in h = a. 
The behavior of the system is not mirrored. 

No matter how we place the steps, all of them inevitably cause oscillation. This will be explained after we study all 
four kinds of external magnetic field. But now we can answer the question: what will happen to the entanglement of 
spins after the constant magnetic field is turned on, and can we benefit from that, say using the step magnetic field as 
an entanglement switch? The answer is negative, because the entanglements oscillate fast and relatively big among 
values. Simple step magnetic field does not provide a way to control or tuning the entanglement in this spin system. 
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FIG. 5: (Color online) Dynamics of C(l,4) in the 7-site system when the step magnetic field is changed from a = 0.5 to fe = 1.5 (before 
the "maximum point" h = 2.61), from a = 3.5 to 6 = 4.5 (after) and from a = 0.5 to fe = 4.5 (big step cross the "maximum point"). 
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FIG. 6: (Color online) Dynamics of C(l,2) in the 7-site system when the step magnetic field is changed from a = 0.5 to fe = 1.5 (before 
the "maximum point" h = 2.46), from a = 3.5 to 6 = 4.5 (after) and from a = 0.5 to fe = 4.5 (big step cross the "maximum point"). 
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B. Exponential magnetic field 



The second kind of time-dependent magnetic field we will look at is the exponential one represented by 

It is a more general form than the step function; when w — > cx), the exponential function turns into a step function. 
Fig. [5] highlights the effect of oj on the concurrences C(l,4) and C(l,2). As one can see the concurrence has the 
trend of following the shape of exponential function where it increases suddenly as the higher value is turned on 
and reaches certain asymptotic equilibrium value. As the transition parameter uj increases the concurrence shows an 
oscillation which increases rapidly as we further increase uj where in that case the behavior resembles the step function 
case. It is interesting to see that the concurrence asymptotic equilibrium value, under the exponential magnetic field, 
coincides with the concurrence corresponding to the higher constant magnetic field, h = b, (the blue (red) straight line 
corresponding to C(l,4) (C(l,2)), which means that the two-dimensional Ising spin system shows an ergodic behavior 
in contrary to the one-dimensional system. In all cases the edge concurrence C(l,2) is higher than the central one 
C(l,4) as expected as the latter shares the entanglement with higher number of other spin pairs. Various a and b 
combinations are tested just as what we did for the step function with similar results to that in Fig. [5] Of course, 
changing the value of the final magnetic field lead to a change of the equilibrium value of the concurrence following it 
up and down. The concurrences C(l,5) and C(l,7) show very similar behavior as C(l,2) and C(l,4) as shown in Fig. 
[5] but with smaller magnitude. 
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FIG. 8: (Color online) Dynamics of the concurrences C(l,2) and C{1,4) in applied exponential magnetic fields of various frequencies 
UJ = 0.1, 0.5 & 1, with a = 1, 6 = 2. 



In the sense of tuning entanglement, exponentially changed magnetic field can be used to vary the entanglement 
from one value to the other smoothly as long as long as its transition rate is slow enough. For instance w — 0.1, that 
is ten percent of the interchange coupling J in energy scale, is a good choice to accomplish this task as we will explain 
latter. 
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FIG. 9: (Color online) Dynamics of the concurrences C(l,5) and C{1,7) in applied exponential magnetic fields of various frequencies 
uj = 0.1, 0.5 & 1, with a = 1, 6 = 2. 



C. Hyperbolic magnetic field 

As we have seen in the previous two subsections, applying a step function or a rapid exponential function may disturb 
the system and lead to a strongly oscillating concurrence. Now let us apply another form of external magnetic field, 
namely hyperbolic which provides a smoother interaction with the system and reduces disturbances. The hyperbolic 
field is represented by 

fa t < to 

= I (6_a) [tanh(L.t) + 1] + a t > to ^^^^ 
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FIG. 10: D(Color online) Dynamics of the concurrences C(l,2) and C(l,4) in applied hyperbolic magnetic fields of various frequencies 
oj = 0.1, 0.5 & 1, with a = 1, 6 = 2. 



Figure. [TOl explores the behavior of the concurrence C(l,4) and C(l,2) of the system under the hyperbolic field at 
different transition constant values uj. Comparing the behavior of the concurrence in this case with the exponential 
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FIG. 11: (Color online) Dynamics of the concurrences C(l,5) and C(l,7) in applied hyperbolic magnetic fields of various frequencies 
uj = 0.1, 0.5 & 1, with a = l,b = 2. 



case at the same frequencies, it is clear the similarity but with much higher controllability for the hyperbolic field 
with sharper asymptotic concurrence value. The system confirms its ergodic behavior as can be seen. Again the next 
nearest concurrences C(l,5) and C(l,7) show the same behavior as the nearest neighbor concurrences as shown in 
Figure. [TTl 



D. Periodic magnetic fields 

In this section we test the dynamics of the system under a different form of external magnetic field, namely periodic. 
It is represented by the sinusoidal function form 

= • f (28) 

^ ' [ a — a sm(u}t + (p) t > to 

where (j) represents an initial phase of the function, which determines the initial value of the applied magnetic field. 
For (/) = we obtain a sm{ujt) function, which is studied in Figs. [T2l and fT3l while for (/i = 7r/2 we get a cos{ujt) which 
is explored in Figs. [HI and[T5l 



The influential factors in the sinusoidal fleld are amplitude a and angular frequency w. Figure 12(b) is a good start 



to analyze them. When ui — 0.1 and a — 1, both are small, concurrence varies up and down in the same frequency 



as the field, but little dents and big catches come out as both w and a increase as shown in Figs. 12(c) and 12(d) 
As the frequency of the field increases it becomes too fast for the system to follow causing more imperfection in the 
concurrence oscillation. On the other hand, as we have seen in the previous magnetic field forms, larger magnetic 
fields don't necessarily bring larger concurrence, so forth dents appear as we increase a. Even larger magnetic fields 



totally break the pattern as can be seen in Fig. 12(d) It is interesting to see that the larger amplitudes are not so 
disturbing to the nearest concurrences as they are for the next nearest neighbors as can be concluded from Fig. 1131 

The critical effect of the initial phase, which determines the initial value of the field, can be seen in Figs. [14] and 
[in] where the phase was chosen to be 7r/2 resulting a cos{ujt) magnetic field. As one can notice the magnetic field is 
initially zero and the concurrence is closely following the magnetic filed with much sharper and much less distorted 
profile than the sin case. In fact, testing a middle value for the phase, between and 7r/2, showed a concurrence 
with a middle profile between the sin and cos cases (Fig. [16]). Therefore, the smaller the initial value of the external 
magnetic field, better if zero, the less distorted is the concurrence profile. 



IV. DYNAMICS OF THERMAL ENTANGLEMENT 



In this section we intend to take a glance at the properties of entanglement in many-body system at finite temper- 
atures. The states describing a system in thermal equilibrium state at absolute temperature T, are determined by 
the Hamiltonian of the system and the inverse temperature /3 = 1/fcT, where k is Boltzmann constant. The thermal 
density matrix of the system is p = Z~^e~^^ , where Z = tr{e~^^) is the partition function of the system. Consider 
our step system and micro-step system shown in Figs. [T]and[21 where any form of magnetic field can be represented 
by a sequence of step functions. 
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FIG. 13: (Color online) Dynamics of the concurrences C(l,5) and C(l,7) in applied sine magnetic fields of various frequencies and field 
strength ui = 0.1 &i a = 1, lo = 0.5 &; a = 1. 



Assuming that the magnetic field has been divided into a sequence of steps, a^b,c, . . . at times tQ,ti,t2, • ■ ■ , then 
at t — to, the system is in an initial equilibrium state described by, 



H{a) = - afa^-aE^at, (29) 

and 



<»j> 



P^^O) - ^ g_;3£;,(a) • i-^") 

At ti > to, the system evolves under the new magnetic field b such that 

H{b) = -J2 ata^^bJ:,a!, (31) 

<h3> 
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FIG. 14: (Color online) Dynamics of the concurrences C(l,2) and C(l,4) in applied cosine magnetic fields of various frequencies and field 
strength a; = 0.1 k a = I, to = 0.5 & a = 1, w = 0.5 & a = 5. 
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FIG. 15: (Color online) Dynamics of the concurrences C(l,5) and C(l,7) in applied cosine magnetic fields of various frequencies and field 
strength u} = 0.1 &c a = 1, uj = 0.5 & a = 1. 



Similarly at t2 > ^i, we have 



and 



P{t2 



<i,3> 



-iH{c){t2-tt) p^^^yiH(c)(t2-tt) 



(32) 



(33) 



(34) 



Continuing in the same way along this sequence we can obtain the density matrix at any time which can be used 
to evaluate the concurrence as explained earlier. 

One can see from the formulation that the thermal equilibrium (relaxation) only enters at t = 0, where we assume 
that the time scale of the dynamics under study is much smaller than the thermal relaxation time. 
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FIG. 16: (Color online) Dynamics of the concurrences C(l,2) and C(l,4) in applied sinusoidal magnetic fields 
frequencies and field strength u) = 0.1 & a = 1, w = 0.5 & a = 1, w = 0.5 & a = 5. 
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FIG. 17: The time evolution of the concurrence as a function of the temperature kT under an applied step magnetic field, where time t 
is in the unit of and kT in the unit of J. 



Figure [T7l shows how concurrence evolves versus time and temperature under step magnetic field h[t) = a aX t < U)^ 
and h{t) = b at t > t^. Adjusting the step values, lead to similar behaviors. They all oscillate through time which 
is consistent with the results at zero temperature. When the temperature increases, oscillations keep the shape but 
are weakened obviously. As can be seen for either the step a = 1, 6 = 2, or the step a — 2, b ~ 3, the concurrence 
disappears around kT = 1.75. kT = 1.75 means when the energy associated with temperature is 1.75 times of the 
exchange interaction J of the system, concurrence will be "killed" , which shows how fragile the entanglement is in 
this spin system. Figures [TS] and \W\ show a very similar behavior of the concurrence of the system under other forms 
of the external magnetic field namely exponential and hyperbolic respectively. The general profile of the concurrence 
resembles the zero-temperature case specially at low temperatures. Remarkably, as the temperature increases the 
concurrence decays very rabidly close to the zero temperature as shown in Figs. 18 (a) [an d 19(a) before it completely 
vanishes at about the same temperature as the exponential case (Figs. [18(b)] andptb)]). 



To understand the critical temperature, at which the concurrence vanishes, better we draw the graph "concurrence 
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C v.s. temperature kT under certain constant magnetic field strength a" (Fig. lIVp . As the magnetic field gets 
big, the critical temperature gets big too. When a is small, zero temperature concurrence is small, and stays stable 
as temperature increases a bit, but for not long it drops dramatically. When a is close to the "maximum point" 
(Fig. [5]), the zero temperature concurrence is relatively large, and drops all the way down as temperature increases. 
When a passes that point, though the concurrence starts smaller and smaller, it owns the ability to be more and 
more stable over temperature change, and finally become zero at relatively larger temperature. This is telling us that 
sacrificing certain amount of concurrence (choosing a large a) may make the concurrence more robust against thermal 
fluctuations. 

In the attempt to understand the vanishing of concurrence, we found that if only the ground state doublets are 
included (others excluded by force), which are not exactly degenerate due to finite size effect, the entanglement never 
vanishes at finite temperature. However, after the one-particle excitation continuum states are included, with higher 
energy states projected out by force, the entanglement disappears at certain temperatures. This tells us that the 
mixing of the ground state and the low-energy excited states will cause the system lose its entanglement. Physically, 
the population of the excited state n is determined by the ratio {En — Eq)/T, and remains almost zero for T < En — Eq 
and becomes significantly nonzero as T > En — Eq. Therefore in the well gapped regions (very small fields and large 
fields) the entanglement sustains the T = values up to a certain T* and then drops to zero; while in the critical 
region the entanglement starts dropping immediately after T > without a plateau region. Numerically, however, T* 
is, thought related, not simply equal to the excitation gap. This might be so because the size of the system prevents 
it from revealing what only statistically holds true. Another question to be studied is why the concurrence should 
come to zero in a 'sudden' way like undergoing a phase transition, in contrast to an exponential decay as one expects 
with a crossover. 
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a=0.1 



FIG. 20: Concurrence vs temperature at certain constant magnetic field strength a, where a varies from 0.1 to 15 in the unit 
of J. 

V. FERMI'S GOLDEN RULE AND ADIABATIC APPROXIMATION 

For all exponential, hyperbolic and periodic, when the transition constant (frequency) uj is small, entanglement 
tends to follow the change of external magnetic field; when uj gets larger, entanglement gradually loses pace with 
the field. These phenomena can be explained by Fermi's golden rule and the adiabatic approximation. In order for 
the entanglement to follow the change in the external field, one requires that the system not to deviate far from the 
ground state, which is the adiabatic approximation. 

Physically, it is equivalent to the requirement that the characteristic frequency in the external field is much smaller 
than the gap, Ei(t) — Eo{t). This may be demonstrated in the following way. From time-dependent perturbation 
theory, given a system at its ground state at t = 0, the probability amplitude of the nth state at t > is given by 



Cn{t) 



dt{n\H'{t)\0) exp{i{En - Eo)t). 



The transition probability from ground state to the nth state is 

P0n{t) = |c„(i)P « %^|/i(i?n - Eo)\', 



where we have 



H'{t) = h{t)Y^al^h{t)S^ 



(35) 



(36) 



(37) 



and 



h{uj') = / h{t) exp{iu;'t)dt. 



(38) 



This is the Fermi's golden rule, that the system only absorbs perturbations at frequencies that match the excitations 
energies. Both data shown for the exponential form and the harmonic form of external field can be explained using 
this principle. For example, if one has 



h{t) 



1 

1 - e" 



t < 
t>0 ' 



(39) 



Fourier transform gives 



'M2 



\Kco')\ 



(cj'2 +a;2)- 



(40) 
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The transition probability is 



From this formula we can see that if one has cj 'ti Ei — Eq, the transition rate is very small and the system is able to 
follow the change of the magnetic field. 

The validity of the adiabatic approximation may fail in two cases. First, when the external field is changing too 
fast for the system to follow, say, a; 3> En — Eq for certain n, and the system has a significant probability of being 
excited to the nth state thus losing the ground state entanglement. One can see this type of breaking for exponential 
(hyperbolic) form of external field with large w and for step function (which corresponds to a; oo of the exponential 
form) and for the case of harmonic field with high frequency. The second type of adiabatic approximation breaking 
takes place if the strength of the external field is too strong that during its change, the system has to cross its phase 
boundary. Suppose the system is able to follow the field and stays on the ground state up to the critical field, then 
right at the transition the gap between the ground state and a one-particle excitation continuum is closing therefore an 
arbitrarily slow field can significantly send the system to various excited states. Although this picture is in principle 
more relevant in the thermodynamic limit than to the small sized system in the current study, this type of adiabatic 
breaking is still observed in Fig. 12(d)| & |14(d)[ 



VI. EXTENSION OF PROJECTION METHOD TO LARGER SYSTEMS 



At the end of section "Step by step projection" , we mentioned the summation of eq. (IT^ covers all the eigenstates. 
For large systems, such as 19-site system, this step is not realistic. Now if we consider a system at zero temperature 
and at i = one starts to change the field. One can prove that not all excited states need to be included in the 
projection, because most excited states have zero overlap with the ground state by symmetry. The global symmetries 
of the system are: spatial six- fold rotation Ce, spin-flip operation Z2, spatial reflection about the x-axis mx and the 
y-axis my. Group theory tells us that the Hamiltonian can be block diagonalized and each block corresponds to 
an irreducible representation of the symmetry group. If the ground state lies in the i-th block, then since the time 
evolution operator U{t) = T exp{—i H{t)dt) commutes with all symmetry operators, the end state |V'(i)) still lies in 
the same irreducible representation. Let us illustrate this idea by considering a simpler symmetry group constituted 
by only Ce and Z2, then the group has twelve one-dimensional irreducible representations. Each representation 
corresponds to the eigenvalues of the two operators: ((— 1)™, exp(i^)), where m = 1,2 and n — 1,2. ..6. Suppose 
IV'o) belongs the sector denoted by (m, n), and \(t>m' n' [t)) is an eigenstate of H{t) in the (m', n') sector. Then we have 

[C6,;7(i)] =0, (42) 



(i) I [C^6,C/(i)] I V'(O)) = (e™V3_e™V3)(^^,^,|^(i)) =0. (43) 

Therefore the overlap is nonzero only if m! = m. Similarly we can prove n' = n. The simplified symmetry group 
can divide the original Hilbert space into (unequally) twelve parts, and the dynamics is fully captured in only one of 
them; the additional rrix.y symmetries can help further reduce the number of states one need to consider for the zero 
temperature dynamics. 

Although symmetry has helped largely reduce the calculation cost, the remaining problem may still be too extensive 
to solve. It is reasonable to apply an approximation: only states that have energies not much higher than the ground 
state are included in the projection if the external field does not vary too fast. Consider the case when the field 
changes very slowly, hence the adiabatic approximation is valid and one needs to include only the ground state in the 
calculation. Further, if the field changes faster than the excitation gap but not as fast as ranging the whole excitation 
spectrum, it would be enough if one only includes the first few excited states (possibly the first band of excitations) 
to cover the most important time scales. 



VII. CONCLUSION AND FUTURE DIRECTIONS 



We have investigated the dynamics of entanglement in a two-dimensional triangular Ising spin lattice in an external 
time-dependent magnetic field. The spins are coupled to each other through nearest neighbor exchange interaction. 
We studied nearest neighbor and next nearest neighbor concurrences of the system under different time-dependent 
forms of the external magnetic field; step, exponential, hyperbolic and periodic. In contrary to the one dimensional 
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Ising spin system, the two-dimensional system shows an ergodic behavior under the effect of the time dependent 
magnetic fields. The step magnetic field causes great disturbance to the system and leads to rapidly oscillating 
concurrence with amplitude and frequency that depend on the magnetic field step values. The System shows more 
controllability under the effect of the other forms of magnetic fields where the concurrence profile follows the shape of 
the applied magnetic field very closely particularly for small magnetic field strength and small transition constant for 
the exponential and hyperbolic fields and frequency for the periodic field. As the values of these parameters increase 
the concurrence breaks the pattern and rapid oscillations take place. The initial value of the applied periodic field, 
independent of the amplitude, is very critical to the oscillating concurrence profile, smaller initial value yields less 
distorted oscillation. Studying the entanglement at zero and finite temperature revealed that it sustains the same 
profile under the difi'erent magnetic fields as the temperature increases but with reduced magnitude. The eff'ect of the 
temperature is very devastating to the entanglement of the system which decays rapidly as the temperature increases. 
Interestingly, though small value of the magnetic field strength leads to small concurrence, it is found to be more 
robust against thermal fluctuations than the larger field strength. In future, we would like to study a larger size 
two-dimensional spin lattice to examine the effect of the size on the different properties of the system [s^ . Also we 
would like to study a more generalized spin system where the coupling among the spins in the other directions, rather 
than the z direction, is taken into account as well as the effect of a time dependent exchange coupling. 

Furthermore the dynamics of a real system is determined not only by its internal Hamiltonian but also by its 
environment. The rich and varied physics of spin systems make spin baths fundamentally interesting [37|. We plan to 
study the reduced dynamics of the center spin with the rest as environment, also entanglements and decoherence for 
a reduced system with two interacting spins and the rest as the environment, in the 7 sites and 19 sites spin systems. 
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